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O - ABSTRACT 

in ; 

Chemical evolution models for the Galactic disk under an inside-out formation 

O . 

scenario are presented for seven sets of stellar yields, all of them metal dependent. In 
particular, the effects of yields from massive and low-and-intermediate-mass stars on 
the C/O chemical history of the solar neighborhood and on abundance gradients are 
discussed comparing the predictions with abundance ratios from nearby H II regions, 



B-stars, dwarf stars, and the Sun. 

In the solar vicinity, the increase of C/O with metallicity is due to massive stars 
alone. Models with yields by Maeder (1992) or by Portinari, Chiosi, & Bressan 
(1998) reproduce the observed C/O increase with metallicity, while models assuming 
yields by Woosley & Weaver (1995) and Woosley, Langer, & Weaver (1993) do not. 
Models based on yields by Maeder (1992) are in agreement with the C/O gradient, 
while those based on yields by Portinari et al. (1998) or by Woosley & Weaver 
(1995) and Woosley et al. (1993) are not. To match both, the C/O enrichment 
history of the solar neighborhood and the C/O local gradient, models call for: i) 
yields that take into account metal- dependent stellar winds; and ii) for massive stars, 
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a mass- loss rate with a metallicity dependence not as simple as Z - 5 . 

C/O abundance ratios predicted with yields by Marigo, Bressan, & Chiosi (1996, 

1998) are higher than those with yields by Renzini & Voli (1981), which in turn are 

higher than those with yields by the van den Hoek & Groenewegen (1997). 

The observed O/H and surface density distributions are matched using any yield set. 

AY/ AZ values are computed for different combinations of seven sets of stellar yields. 

All models agree with the AY j AZ value derived from M17, but are smaller than 

those derived from other H II regions. 

Key Words: GALAXY: ABUNDANCES - GALAXY: EVOLUTION - 
STARS: EVOLUTION 

1. INTRODUCTION 

In a chemical evolution model that assumes an infall of H and He gas without mass 
loss by galactic winds or radial flows, the predicted C/O abundance ratio depends 
mainly on the stellar yields and on the initial mass function. By fixing the initial 
mass function, this work explores the use of C/O to test stellar yields. 
The stellar yield of an element Xj, being the mass fraction of a star of initial mass m 
converted to X i and ejected to the interstellar medium (ISM), strongly depends on 
the assumptions of the stellar evolution modeling. By the 80', no stellar yield set for 
massive stars (m > 8 M ) of different metallicities had been computed. In this 
decade, Maeder (1992, M92) published He, C, O, and Z stellar yields for two initial 
metallicities (Z = 0.001 and 0.02) for objects with initial mass in the 9-120 M 
range. Later on, Woosley & Weaver (1995, WW) calculated stellar yields for 80 
elements and isotopes for stars with initial masses between 11 and 40 M and with 
five initial metallicities (Z/Z e = 0, 10~ 4 , 0.01, 0.1, and 1), without initial abundance 
ratios scaled to solar. More recently, Portinari, Chiosi, & Bressan (1998, PCB) 
published stellar yields of 17 elemental species for a wide range of metallicities 
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(Z = 0.0004, 0.004, 0.008, 0.02, and 0.05) and masses (6-120 M ). 

For low and intermediate mass stars (LIMS, 0.8 < m/M < 8), stellar yields have 
been calculated by Renzini & Voli (1981, RV), van den Hoek & Groenewegen (1997, 
HG), and Marigo, Bressan, & Chiosi (1996, 1998, MBC). While HG produce a grid 
of stellar yields for five metallicities (Z =0.001, 0.004, 0.008, 0.02, and 0.04), RV and 
MBC sampled yields for only a pair of metallicities [Z =0.004, 0.02 and Z =0.008, 
0.02 respectively). A combination of PCB and MBC provides the only complete grid 
of stellar yields computed from the same evolutionary tracks (Padova group). 

Carigi (1994) has shown that the evolution of [C/O] with [O/H] in the solar vicinity 
can be explained by M92 yields due to the dependence of the C and O yields on 
initial Z. Prantzos, Vangioni-Flam, & Chauveau (1994) concluded that the 
metallicity dependent M92 yields are also able to reproduce the growth of [C/O] in 
the solar neighborhood. Extending the modeling beyond the solar vicinity, this work 
follows the [C/O] abundance and evolution with chemical evolution models for the 
galactic disk, with the aim to study the behavior of differences of sets of 
metal-dependent yields from both, massive stars and LIMS. 

In general, chemical evolution models of the solar vicinity and the Galactic disk 
assume a given yield set, either M92 yields (e.g. Giovagnoli & Tosi 1995, Prantzos & 
Aubert 1995, Carigi 1996), WW yields (e.g. Timmes, Woosley, & Weaver 1995, 
Chiappini, Matteucci & Gratton 1997, Allen, Carigi, & Peimbert 1998, Prantzos & 
Silk 1998), or Padova yields (e.g. Portinari et al. 1998, Tantalo et al. 1998). Since 
these yield sets have substantial differences, it is not straightforward to separate the 
effects of yields from other model assumptions when intercomparing the different 
studies. In this work, consistent chemical evolution models are build to find out 
which yield sets reproduce better the ISM abundances. The predicted abundances 
from the different sets of stellar yields are compared for the four elements in common 
among all sets (H, He, C, and O). 
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In this paper, all chemical abundances are presented by number, with the exception 
of AY/ AO, and AY/AZ, that are given by mass. 

The paper has been organized as follows: The observational constraints on the 
models are presented in §2. §3 describes the assumptions adopted in the chemical 
evolution models. The model predictions are then shown and briefly discussed in in 
§4. The global discussion and conclusions are presented in §5 and §6, respectively. 

2. OBSERVATIONAL CONSTRAINTS 

In this work, the compilation by Peimbert (1999) and the recent data of Esteban et 
al. (1998, 1999a, 1999b) are used as abundance constraints for the models. He/H, 
C/H, and O/H values, for the Galactic Hn regions M17, M8, and Orion are taken 
from Peimbert's Table 1, but corrected for dust depletion, adding 0.08 dex and 0.10 
dex to the O and C abundances respectively, as suggested by Esteban et al. (1998). 
In the modeling, I attempt to fit the abundance gradients to be within the average 
values by Peimbert (1999) and those obtained from the observations by Esteban et al. 
(1998, 1999a, 1999b), considering temperature fluctuations (t 2 > 0.00). Average O/H 
gradients computed by Peimbert (1999) are based on the gradients from Esteban et 
al. (1998, 1999a, 1999b), Shaver et al. (1983), and Deharveng et al. (1999). Since in 
the literature there are C/O values based on recombination lines only for M17, M8, 
and Orion, the average C/H gradient is obtained from Esteban et al. (1998, 1999a, 
1999b) and Peimbert et al. (1992). The relative enrichment AY/ AO and AY/AZ 
values for M17, M8, and Orion are taken from Peimbert (1999, t 2 > 0.00), which 
have already been corrected for dust (also as suggested by Esteban et al. 1998). 
O/H abundances and gradient from B-stars are in agreement with those from Hn 
regions (Gummersbach et al. 1998, Smartt & Rolleston 1997), but the C/H values 
are lower by at least 0.3 dex, and the C/H gradient is twice as flat (Gummersbach et 
al. 1998). The gradients from B-stars shown in Table 3 were computed using the 
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data of Gummersbach et al. (1998) and Smartt & Rolleston (1997) for the 
galactocentric range considered here (4 < r < 10 kpc). The innermost B-star of both 
data sets was eliminated from the fit, given their uncertain abundances as discussed 
by the respective authors. The He/H data by Gummersbach et al. (1998) are not 
considered in this work given that the large dispersion in abundances 
(-1.1 < log(He/H) < -0.5, corresponding to 14 < AY/AO< 150) does not provide 
a tight enough constraint for the chemical evolution modeling. The large differences 
in He abundance observed among B-stars may come from a contamination in the 
outer stellar layers by helium produced by the stars, combined with errors and other 
uncertainties in the abundance determinations. Although B-stars and Hll regions 
show a different slope of the C/O gradient, both of them are negative. In this study, 
the presence (not the slope magnitude) of a negative C/O gradient in the galactic 
disk is used as another observational constraint. 

The model C/O history for the solar vicinity is constrained with C/O abundances 
from dwarf stars located 1 kpc around the Sun (Gustafsson et al. 1999) and the solar 
value from Grevesse & Sauval (1998). Since C/H values from B-stars are lower than 
those from Hll regions, two different present-day C/O values are presented. The 
increase of C/O with metallicity in the solar vicinity provides yet another 
observational constraint. 

Planetary nebulae abundances are not be used here as constraints because C is 
produced by their progenitors. The models by Renzini & Voli (1981), van den Hoek 
& Groenewegen (1997), and Marigo et al. (1996, 1998) predict C enrichment in the 
envelope during the evolution of stars with masses lower than ~ 8 M Q . Moreover, 
Peimbert, Torres- Peimbert, & Luridiana (1995) and Peimbert, Luridiana, & 
Torres-Peimbert (1995) found C/H values higher than the solar value, by at least 0.1 
dex, for the vast majority of planetary nebulae in their sample. 

Another important observational constraint is gas consumption, measured as the 
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ratio of gas to total surface mass densities, o gas jo tot . In this work the observed a gas 
distribution, as compiled by Matteucci & Chiappini (1999), is used together with an 
exponential a tot distribution with a 3 kpc scale-length and a cr tot {r- & ) = 45 M Q pc~ 2 
amplitude (Kuijken & Gilmore 1991). A galactocentric distant for the Sun, r Q , of 8 
kpc is adopted here. 

Summarizing, the just-described observational constraints are used in this study as 
follows: a) All models (for each and every yield set) are build to exactly reproduce: 
i) the observed gas fraction distribution of the galaxy, o gas j 'atot] and ii) the observed 
O/H galactic gradient, b) In order to study the differences among available yields 
sets, the model predictions are then compared to: i) the observed rise of C/O with 
metallicity (or equivalently with time) in the solar neighborhood; and ii) the 
observed decrease of the C/O abundance with galactocentric distance. 

3. CHEMICAL EVOLUTION MODELING 

The scope of this work is to explore models assuming the different sets of 
metal-dependent stellar yields available in the literature. 

All models are built to reproduce two observational constraints: the O/H abundance 
gradient from Hll regions and B-stars, and the a gas distribution from 4 to 10 kpc. 
Models are very similar to the infall model of Allen, Carigi, & Peimbert (1998), but 
adopting somewhat different assumptions: 

a) The star formation rate is set proportional to a power of a gas and a tot : 
SFR(r,t) = v <Jg as (r,t) af 7 1 (r, t), where v and x are constant in time and space, and 
fixed in such away that the observational constraints are reached after 13 Gyr, the 
age of the model. 

b) Three sets of metal-dependent stellar yields from massive stars are used: i) Geneva 
yields: M92 for 9 < m/M < 85 (high mass-loss rate); ii) Santa Cruz yields: WW for 
11 < m/M Q < 40 (models "B" for 30, 35 and 40 M ), and WLW for m = 60 M and 
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m = 85 M ; Since Z max = Z Q for the Geneva and Santa Cruz groups, then yields 
(Z > Z Q ) = yields(Z ). iii) PCB yields: Portinari et al. (1998) for 9 < m/M < 85. 
In all cases, effects due to black hole formation are not considered. 

c) Three sets of metal-dependent stellar yields for LIMS are used: i) RV yields: 
Renzini & Voli (1981) for 1 < m/M < 3 (a = 0.0, -q = 0.3, case A), for 

3 < m/M < 8 (a — 1.5, rj = 0.3, case A), ii) HG yields: van den Hoek & 
Groenewegen (1997) for 0.8 < m/M < 8 (jjagb = 0A,M HBB = 0.8). iii) MBCP 
yields: Marigo et al. (1996) for 0.8 < m/M < 3; Marigo et al. (1998) for 4, 4.5, and 
5 M ; Portinari et al. (1998) for 6 and 7 M . 

d) The fraction of binary systems predecessors of SNIa is taken as A = 0.07. 
Contrary to Carigi (1994), SNIb produced by binary systems are now not considered. 

Carigi (1996) and Allen, Carigi, & Peimbert (1998) successfully modeled galactic 
abundance gradients considering the IMF by Kroupa, Tout, & Gilmore (1993). 
Therefore, in this work the same IMF is assumed for a 0.01 to 85 M mass range. 
Moreover, this IMF is maintained constant in time and space, based on the 
conclusion of Chiappini, Matteucci, & Padoan (2000) that a constant IMF is still the 
best assumption to explain the observational constraints in the Milky Way. 

In the present exercise, to avoid spurious mass extrapolation effects for a fair 
comparison of the predictions from the different yield sets, the stellar mass range was 
limited so as to be lower than 85 M (the maximum mass common to the three sets 
from massive stars). 

Since the code does not use the instantaneous recycling approximation, it is 
necessary to estimate, at each time step, the stellar yields of all dying stars at the 
metallicities at their birth times. Therefore, the stellar yields from the different 
groups have been interpolated in both mass and metallicity, taking in consideration 
their different samplings. In particular, given the number of initial metallicities 
considered by stellar evolution models, and the yields Z-dependencies, I have 
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assumed: a) a linear interpolation with mass for yields of massive stars and LIMS; b) 
a linear interpolation with metallicity for the Santa Cruz, PCB, RV, HG and MBCP 
yields; and c) three kinds of Z interpolations (linear, as Z a5 , and Z 2 ) for the Geneva 
yields. It is important to remark that the HG yields vary almost linearly with Z, but 
the Z-behavior of yields of massive stars is not well determined: the PCB yields 
dependency on Z is not simple and the Santa Cruz yields are almost independent of 
Z. 

Input parameters for the set of models are summarized in Table 1. Column (2) lists 
x, the SFR power on a gas , and column (3) the efficiency factor, v, of each model. 
The slope of the final abundance gradients is basically governed by x; while v and x 
together determine the O abundance. Models are divided in three groups (column 
4), assuming M92, WW&WLW, or PCB yields. M92 and WW&WLW yields have 
been complemented with RV, HG, and MBCP yields. Three subgroups of models 
have been computed according to the type of interpolation used (Z n , n =0.5, 1, or 2) 
between the two metallicities considered by Maeder (1992). The table lists just one 
model for the PCB yields (C-model), since the Padova yields (PCB+MBCP) form a 
complete and consistent set. 

Since LIMS do not produce O, each set of yields from massive stars requires different 
values of x and v to match the same observational constraints. Since the M92-0 
yields for m > 25 M Q decrease with metallicity, O gradients saturate quicker, calling 
for higher x values. On the other hand, WW-0 yields are almost independent of 
metallicity, and WLW-0 yields are lower than M92-0 yields at low metallicity, 
forcing higher v values for the WW&WLW modeling to reach the observed O/H of 
Hll regions. Since PCB-0 yields are lower than M92- yields at low Z, it is required 
that fc-modci ~ ^M-modeis when the interpolation is weighted to high metallicities. 

Since the yields for supermetallic stars are similar to those by WW, then 

^C-model ~ ^W-models- 
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TABLE 1 



MODEL INPUT PARAMETERS 


Model 


x a 


v (Gyr- 1 (M pc" 2 ) 2 * 1 "*)) a 


Massive Stars + LIMS Yields b 


n c 


MXO 


1.22 


0.048 


M92 + X 


0.5 


MX1 


1.22 


0.040 


M92 + X 


1.0 


MX2 


1.22 


0.033 


M92 + X 


2.0 


WX 


1.13 


0.075 


WW&WLW + X 




C 


1.15 


0.045 


PCB + MBC 





a SFR = v a* gas af- t l 
b X = RV, HG, MBCP 
c Interpolation type (Z n ) 
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Since infall models of the present kind produce abundance gradients that flatten out 
with time (Carigi 1996), the M- models require higher x exponents to reproduce the 
O/H gradient at a fixed age (13 Gyr). 

4. MODEL RESULTS 

All models match well several observational constraints at the solar vicinity, e.g. 
age-metallicity relation, G-dwarf distribution, [C/Fe] and [O/Fe] vs [Fe/H] evolution. 
For a SFR= v a x gas af^ 1 law, the present models adopt x in the range 1.13 and 1.22. 
Predictions for the solar vicinity do not change significantly when the x exponent 
varies from 1 to 2 (Carigi 1996). This work does not then discuss how well the 
observed relations are matched, since the predicted relations are very similar to those 
of the best model of Carigi (1996), where her Fig 4 already summarizes the 
comparison with the observational data. 

4.1. C/O evolution of the solar vicinity assuming different yields 

Figure 1 shows the C/O evolution in the solar vicinity as predicted by models that 
consider yields by Maeder (1992). Since the dependence of these yields for massive 
stars with Z is not well known, the influence of the assumed type of interpolation on 
on the C/O ratio is explored in Figure la. The figure shows that, between 1.5 and 9 
Gyr, the C/O rise in time slows down as the interpolation exponent n increases (Z n ); 
out of this lapse, the increase is basically the same for all reasonable ns. It is then 
important to compute yields for intermediate metallicities to reduce these 
interpolation ambiguities. Besides interpolations effects, Figure lb also explores the 
influence of LIMS on the C/O evolution, since these stars are important carbon 
producers. After 1 Gyr or so, RV predict less carbon production than MBC but 
more than HG. 

The best models from the M92 yields among the different combinations of LIMS 
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Fig. 1. The C/O evolution of the solar vicinity. Model predictions are presented as follows: (a) 
Considering three different interpolations (Z - 5 , Z 1 , and Z 2 ) of the massive-stars yields by Maeder 
(1992) and LIMS yields by Renzini & Voli (1981). Long- dash- dot lines: MRVO-model, continuous 
lines: MRVl-modcl, short- dash- dot lines: MRV2-model. (b) Assuming the Z 1 interpolation for yields 
by Maeder (1992) and three set of LIMS yields by: Marigo et al. (1996, 1998), Renzini & Voli (1981), 
or van den Hoek & Groenewegen (1997). Long-dash lines: MMBCl-model, continuous lines: MRV1- 
model, short-dash lines: MHGl-modcl. Observational data are as follows: filled circle: computed 
value at r — r Q from radial gradients by Peimbert (1999); filled triangle: average value for the five 
B-stars at r = r Q + 0.4 kpc from data by Gummersbach et al. (1998); filled squares: dwarf stars at 
r = r Q ± 1 kpc from Gustafsson et al. (1999); 0: solar value from Grevesse & Sauval (1998). The 
ages of the dwarf stars were scaled to the age of the models. 
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yields and interpolation types assumed (M-models = MRV1, MHG1, MMBC2) are 
shown in Figure 2. Results from models assuming the Santa Cruz and Padova yield 
sets are also presented in similar diagrams (Figure 3 and Figure 4 respectively). 
These results, from the modeling assuming the different yield sets, are discussed in 
the next two subsections, looking first at the C/O evolution in the solar vicinity, and 
later on at the local chemical distributions along the galactic disk. 

In particular, the (h) panel of each figure presents the predicted C/O chemical 
enrichment history of the solar neighborhood, as compared to the following 
observational data: i) the observed abundance in the Sun and in the neighboring 
disk dwarf stars, ii) a C/O value at r & derived from the observed gradient for Hn 
regions, and iii) a C/O average value for the nearest B-stars from Gummersbach et 
al. (1998). From these panels, one can see that the C/O abundance always increases 
with time in the M-models and the C-model, in agreement with most of the 
observations (dwarf stars, Sun, and Hn regions), but the results of the W- models fail 
to predict the observed recent rise in C/O. 

The discrepancy among the C/O predictions with different yields can be understood. 
The M92 and PCB-0 yields increase with mass and decrease with metallicity, while 
the C yields increase with both mass and metallicity for Z < Z Q . At the onset of the 
evolution, the WLW yields from very massive stars (m > 40 M ) enrich the ISM 
with ever more oxygen through O-yields that increase with mass. Once these stars 
die after a few million years, the C/O of the W-models behaves like in the M-models 
because, for m < 40 M , both sets of yields are similar at low metallicities. Since the 
WW yields do not depend on metallicity, the W-models reach (at ~ 3.5 Gyr) a 
plateau in C/O that decreases slowly with time. 

The C-model predicts two plateaus in the C/O evolution of the solar neighborhood: 
a first one at ~ 5 Gyr and another after ~ 11 Gyr. The later one is due to the 
contribution of supermetallic stars. These stars have extreme winds that occur 




Fig. 2. a) Present-day radial distribution of gas surface mass density, the area enclosed by dotted 
lines indicates the data collected by Matteucci & Chiappini (1999); (b)-(g) observed and predicted 
abundance ratios. Predictions from the best models considering yields by Maeder (1992): long-dashed 
lines: MMBC2-model, continuous lines: MRVl-model, short- dash lines: MHGl-model. Observa- 
tional data: filled circles: gas and dust in Hn regions from Peimbert (1999) and Esteban et al. 
(1998), respectively; filled triangles: B-stars from Gummersbach et al. (1998); and crosses: B-stars 
from Smartt & Rolleston (1997). The observed values correspond to M17, M8, and Orion at adopted 
galactocentric distances of 5.9, 6.5, and 8.4 kpc, respectively. Error bars at the left and the right 
represent the typical errors from data by Gummersbach et al. (1998) and Smartt & Rolleston (1997), 
respectively, (h) The C/O evolution of the solar vicinity, observational data as in Fig 1. 
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Fig. 3. Same as previous figure but with the predictions from models assuming yields 
by Woosley & Weaver (1995) and Woosley et al. (1993) for massive stars and Marigo 
et al. (1996, 1998), Renzini & Voli (1981), or van den Hoek & Groenewegen (1997) 
for low and intermediate mass stars. Long-dash lines: WMBC-model, long- short- dash 
lines: WRV-model, short-dash lines: WHG-model. 
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Fig. 4. Predictions from models assuming yields by Marigo et al. (1996, 1998) for 
low and intermediate mass stars in combination with yields for massive stars by and 
Woosley & Weaver (1995) and Woosley et al. (1993), Portinari et al. (1998), or Maeder 
(1992). Short-dash lines: WMBC- model, continuous lines: C-model, long-dash lines: 
MMBC2-model. Observational data points and details as in Fig. 2. 
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before C and are synthesized, so their winds are He rich, while C and are 
produced during the SN explosion. These characteristics of the yields from 
supermetallic stars mimics the behavior predicted by WW. 

Table 2 summarizes the present-day predictions for the solar neighborhood from all 
models. M-models reproduce the observed C/H and C/O values from Hll regions. 
WRV and WHG models predict C/H and C/O values in agreement with B-stars, 
while the WMBC model results agree with the C/H and C/O from Hll regions. 
Inspecting Figure 1 and Table 2 it can be noticed that the average C/O for B-stars is 
lower than for Hll regions and for dwarf stars, mainly due to the low C/H values 
observed in B-stars. 



4.2. Chemical radial distributions in the local disk assuming different yield sets 

Predictions for the Galactic disk are also summarized in the (a) to (g) panels of 
Figures 2-4, together with Table 3. From these figures and table it can be seen that: 

a) The C/H gradients predicted by the best models with M92 yields are similar to 
the observed ones from H II regions. The gradients from the W-models (considering 
WW&WLW yields) or the C-model are flatter by about 0.04 dex kpc -1 and are close 
to the gradient traced by B-stars; 

b) C/H abundance ratios from the M, WMBC, and C models agree with the Hll 
regions observations, but those from WRV and WHG models are lower (by at least 
0.1 dex) in closer agreement with the C/H of B-stars; 

c) M-models are the only ones that reach the high C/H and C/O values of Ml 7; 

d) The negative C/O gradients predicted by the M-models are in very good 
agreement with the observations. The gradients predicted by W-models or C-model, 
being positive, are a poor match to the observed gradient; 

e) The C/O values from W-models agree with those determined from B-stars, while 
the C/O values from M and C models agree with the C/O values determined from 
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TABLE 2 

CURRENT PREDICTED AND OBSERVED ABUNDANCE RATIOS FOR THE 
SOLAR VICINITY 



Model 


12 + loe-fC/H) 


12+loefO/H) 




MRVO 


8.74 


8.74 


+0.00 


MRV1 


8.63 


8.73 


-0.10 


MRV2 


8.48 


8.74 


-0.25 


MHGO 


8.66 


8.75 


-0.10 


MHG1 


8.54 


8.75 


-0.21 


MHG2 


8.38 


8.75 


-0.37 


MMBCO 


8 81 


8 75 


-1-0 06 

1 \J . \J\J 


MMBC1 


8.72 


8.74 


-0.02 


MMBC2 


8.60 


8.75 


-0.15 


WRV 


8.32 


8.71 


-0.39 


WHG 


8.26 


8.71 


-0.45 


WMBC 


8.49 


8.73 


-0.23 


C 


8.62 


8.73 


-0.12 


Obs. M17, M8, Orion a 


8.54 ± 0.12 


8.74 ± 0.07 


-0.20 ± 0.08 


Obs. Hn regions Average a 


8.53 ± 0.12 


8.74 ± 0.07 


-0.21 ± 0.08 


Obs. B-stars b 


8.16 ± 0.12 


8.62 ± 0.08 


-0.46 ± 0.13 


Obs. B-stars c 




8.84 ± 0.10 





a Computed from data by Peimbert (1999) at r 

6 Average values for the five B-stars at r = r + 0.4 kpc from data by Gummersbach 
et al. (1998) 

c Average values at r = r ± 1 kpc from data by Smartt & Rolleston (1997) 
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H II regions; 

f) M-models and the WMBC-model reproduce the observed He/H abundances in Hll 
regions, with the former yielding the best agreement; 

g) RV yields predict lower C and He than MBC but higher than HG; 

h) M-models predict strong negative AY/ AO gradients while W-models and 
C- model yield almost no gradient; 

i) Predicted AY/ AZ gradients are almost flat for all models; 

j) With the exception of the WMBC model, all models predict helium to 
heavy-elements abundance ratios in the 1.6 < AY/ AZ< 1.0 range, values that are 
lower than those observed in M8 and Orion. 

The C/O gradients from the W-models and the C- model are almost flat because 
WW yields are almost independent of Z and the very metal-rich stars eject C and O 
in SN explosions and not in winds, both facts imply a saturation in C/O. 
The He/H gradients predicted by the M-models are steeper than by the C-model 
because the very-massive and super-metallic stars eject a smaller amount of He than 
solar metallicity stars (PCB) 

4.3. C/O vs O/H 

In Figure 5 the predicted evolutions of C/O versus O/H for the solar neighborhood 
are presented. In addition to Galactic data, C/O and O/H values for extragalactic 
H II regions are shown. Objects in the halo and in the disk of the solar neighborhood 
(cylinder of radio 1 kpc centered in the Sun), and in external galaxies call for an 
increase of C/O with O/H enrichment. Only models that reproduce this increase are 
shown in the figure, so the W-models are excluded. The agreement for log(0/H) < 
-3.8 is not so good, with the exception of MHG1. Again, here one can notice that C 
yields for intermediate mass stars given by HG are lower than those obtained by RV 
or MBCP. 
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TABLE 3 

PREDICTED AND OBSERVED PRESENT-DAY RADIAL GRADIENTS 
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-0.038 


-0.053 


+0.014 


C 


-0.056 


-0.054 


-0.002 


Obs. H II regions Average c 


-0.106 


-0.055 


-0.051 


Obs. B-stars d 


-0.059 


-0.038 


-0.026 


Obs. B-stars e 




-0.017 





a Between 6 and 8 kpc 
b in dex kpc -1 
c Peimbert (1999) 

d Computed from data by Gummersbach et al. (1998) (r — 2.5 kpc < r < r© + 1.8 
kpc) 

e Computed from data by Smartt & Rolleston (1997) (r Q — 2.5 kpc < r < r + 1.9 
kpc) 
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The log(C/0)-log (O/H) relation predicted by the MHGl-model of this work (Fig 5) 
is not surprisingly very similar to the one obtained by Henry, Edmunds, & Koppen 
(2000), since their best models, either analytical or numerical, are constructed also 
adopting the yields of Maeder (1992) and van den Hoek & Groenewegen (1997). 
A powerful tool for comparing chemical evolution models with interstellar medium 
abundances is given by the relation 

log(C/0) = a log(0/H) + b. (1) 

Table 4 presents the a values predicted by the models, considering different 
concentric galactic rings: i) a very local disk, between r = 6 and r = 8 kpc, and ii) 
the complete area studied in this paper. Moreover, this table shows the a values 
derived from four sets of observations, two galactic and two extragalactic ones; for 
B-stars it was assumed that the C and O abundances correspond to those of the 
interstellar medium. The observed a values indicate that the negative C/O local 
gradient may be extended to the whole Galactic disk. As can be seen from the table, 
the observed a values are in good agreement with the M-models and in disagreement 
with the W-models and the C-model, reiterating the results presented in Figure 5. 
According to the discussion of Figure 5, the best models are: MRV1, MHG1, and 
MMBC2. 

5. DISCUSSION 

The determinations of C/O from Hll regions in spiral and irregular galaxies 
(Garnett et al. 1999), Hll regions in our Galaxy (Peimbert, Torres-Peimbert & Ruiz 
1992, Esteban et al. 1998, 1999a, 1999b), and dwarf stars in the solar vicinity 
(Gustafsson et al. 1999); all converge to the conclusion that C/O must increase with 
metallicity. M92 yields can explain this behavior, Padova yields only partially, while 
WW&WLW yields not at all. Other authors (Ferrini, Matteucci, Prantzos, Timmes, 
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Fig. 5. Log(C/0) - log(0/H) relation for r = rQ. The predictions of selected models (those that reproduce this 
relation) are presented as follows: dotted line: C-model, long-dashed line: MMBC2-modcl, continuous line: MRV1- 
model, short-dash line: MHGl-model. Observational data for the solar vicinity (filled symbols) arc as follows: circle: 
Orion from Peimbert (1999); hexagons: field halo dwarfs from Tomkin et al. (1992); squares: dwarf stars from 
Gustafsson et al. (1999), as in Fig 2; triangle: average B-star computed from the Gummersbach et al. (1998) data (as 
in Fig. 2); 0: solar value from Grevesse & Sauval (1998). Error bars at the upper-left and the lower-right are typical 
errors from data by Gustafsson et al. (1999) and Tomkin et al. (1992), respectively. Log (O/H) for halo dwarfs is 
the average of [O/H] values from Tomkin ct al. (1992). Observational data for external galaxies (open symbols) are 
as follows: circles: I Zw 18NW and I Zw 18SE, values adopted from Figure 4 of Garnett et al. (1999); hexagons: 
I Zw 18NW and I Zw 18SE, abundance ratios taken from Table 5 of Izotov & Thuan (1999); abundance data for 
spiral galaxies (squares) and dwarf irregular galaxies and Magellanic Clouds (triangles) are obtained from Garnett et 
al. (1999). Plotted C/O values for spiral galaxies are the average of the abundance ratios for two different reddening 
laws. 
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TABLE 4 

PRESENT-DAY RISE ACROSS THE GALACTIC DISK OF C/O WITH O/H 







Predictions 
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0.56 
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0.98 




0.79 


MRV2 




1.41 




1.03 


MHGO 




0.67 




0.56 
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0.81 
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1.27 




1.01 
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0.31 




0.31 
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0.61 




0.50 
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0.67 


WRV 




-0.05 




-0.05 
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0.07 




0.07 
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-0.27 




-0.28 


C 




0.04 




0.06 






Observations" 






Object 


Ar (kpc) 




CtAr 




H II regions Average b 


2.5 




0.93 ± 0.60 




B- stars c 


12.0 




1.69 ± 2.34 




NGC 2403 d 


5.3 




0.50 ± 0.43 




M101 d 


15.1 




0.98 ± 0.29 





a The table presents a values derived from log(C/0) = a log(0/H) + b, and Ar is in kpc 
b Computed from data by Peimbert (1999) 

c Computed from data by Gummersbach et al. (1998) and Hibbins ct al. (1998) 
d Computed from data by Garnett et al. (1999) 
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and Tosi, see Tosi 1996 for references) using different chemical evolution models, 
have predicted flat C/O gradients (Tosi 1996) because all of them have adopted 
either yields at a fix Z or the WW yields, which are almost independent of Z. 

Carigi (1994) and Prantzos et al. (1994) had already concluded that the evolution of 
[C/O] with [O/H] or [Fe/H], respectively, in the solar vicinity can be explained by 
the metal dependent yields of Maeder (1992). The present work have extended these 
previous ones to the Galactic disk using the recent C determinations in H II regions 
and B-stars as observational constraints. Moreover, seven sets of yields have been 
considered here, composing the complete sample of stellar yields dependent on the 
initial stellar metallicity available in the literature. 

A main difference among these sets of yields resides on the stellar-wind assumptions. 
WW do not consider stellar winds at any stage, and M92 and PCB assume a 
mass-loss rate proportional to Z - 5 and m 2 5 during the post-main-sequence phases. 

C and O yields from M92, PCB, and WW are similar for m < 25 M , where stellar 
winds are not so relevant. In M92 and PCB, most of the net C is produced before 
the SN stage, and since the wind strength increases with metallicity, the C yield 
increases with Z. On the other hand, O is mainly ejected during the SN stage and, 
since the amount of O expelled by the SN depends on the stellar mass just before the 
SN explosion, the M92 and PCB-0 yield decreases with Z. This fact causes the 
M- models and the C- model to reproduce the rise of the C/O in the solar vicinity. 

The increase of C yields and the decrease of O yields with Z can not be extrapolated 
to Z > Z Q . According to PCB, the C to O yields ratio decreases from Z = Z & to 
Z = 2.5Z Q . Supermetallic and massive stars have early and intense winds, which 
occur before these stars synthesize C and O. Therefore, their C and O yields are 
similar to those of metal-poor stars that mainly eject C and O in the SN explosion. 
This behavior of the C and O yields produce a saturation of the C/O abundance 
ratio at ~ -0.1 dex, which produces flatter gradients for the C-model. Since 
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M-models assume that the supermetallic stars behave like metal-solar stars, the 
C-model would predict a negative C/O gradient if the mass-loss rate for Z = 2.5Z & 
is similar to that for Z = Z Q 

It is important to note that in the W-models the wind contribution to the C and O 
yields was not included, since WLW present information to calculate only the He 
expelled by winds. If winds were considered in these models, the C yield would 
certainly be higher but, given the low weight at the high-mass end of the initial mass 
function, hardly enough to reproduce the present-day C/O abundances. When the 
WLW yields are not taken into account (by reducing the upper mass-limit to 40 M ) 
the predicted C/O decreases at times earlier than 0.4 Gyr, but still remains low and 
basically constant for the rest of the evolution (Carigi & Peimbert 2000). 

The M-models and C-model agree better with the younger objects in the solar 
vicinity than with the older ones. The agreement would improve if metal-poor stars 
eject more oxygen and less carbon than predicted by a Z ' 5 law. The M-models are 
in agreement with data for spiral and irregular galaxies while the W-models are not, 
as discussed by Garnett et al. (1999). Then again, the general agreement would be 
even better if the C and O yields for metal-poor stars were more dependent on 
metallicity than assumed by M92 and PCB. Furthermore, the C/O value of I Zw 18 
determined by Garnett et al. (1995, C/O— —0.60 ± 0.10 dex) would be hard to 
explain (not so the lower value of Izotov & Thuan 1999, 

< log(C/0) >= -0.78 ± 0.03 >). It should be noted that the models presented in 
Figure 5 were made to reproduce the Galactic disk and therefore, the comparison of 
the extragalactic Hll regions and halo objects is only indicative. Specific models for 
each galaxy or for the Galactic halo should be carried out, see for example the 
models by Carigi, Colin, & Peimbert (1999) and Carigi & Peimbert (2000). 

Moreover, the M-models and the C-model agree well with the [C/H], [C/O], and 
[C/Fe] vs [Fe/H] relations observed in dwarf stars of the solar vicinity by Gustafsson 
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et al. (1999), after correcting for the adopted solar abundances (Grevesse & Sauval 
1998 as opposed to Anders & Grevesse 1989). The predicted relations are matched 
at high [Fe/H] and become slightly lower than those observed at lower metallicities; 
again, a behavior that calls for a stronger metal dependency of the C and O yields 
for metal-poor stars than a Z - 5 law. W-models basically match the relation slopes, 
but again predict a lower C abundance than observed (by ~ 0.3 dex in C/H and ~ 
0.2 dex in C/Fe). 

Maeder (1993), Woosley & Weaver (1995) and Portinari et al. (1998) have assumed 
different values for the C 12 (a,7)0 16 reaction rate. It is known that this rate is 
uncertain by about a factor of 2. This uncertainty causes minor differences in the 
behavior of C/O as a function of metallicity predicted by the models, but more 
severely affects the absolute C/O value. This change in the absolute C/O value does 
not alter the conclusions of this paper. In other words, to a very good approximation, 
a change in the C 12 (o;,7)0 16 rate modifies b in equation (1), but not a. 

The C/H values derived from Galactic B-stars are about 0.3 dex lower than those 
from Galactic H II regions and from disk dwarf stars of the solar vicinity. Carbon 
abundances of B-stars are not easy to understand and may not be representative of 
the present-day C in the ISM, because: i) being so young (10~ 3 - 10~ 2 Gyr), B-stars 
would be expected to be richer in C than most dwarf stars, but the observed C/H in 
B-stars in the solar vicinity is very similar to the C of the poorest dwarfs; ii) the O 
and N abundances and gradients from B-stars and from H II regions are very similar 
(surprisingly so for N), so it is difficult to understand why the C/O gradient from 
B-stars is almost flat, while the one from Hn regions is not; iii) a real C/O gradient 
is supported by the fact that other two nearby spiral galaxies, M101 and NGC 2403, 
also show significant C/O gradients, -0.04 and -0.05 dex/kpc respectively (Garnett et 
al. 1999); and iv) C/O from Hn regions, the Sun, and dwarf stars are consistent 
with each other, but not with the C/O from B-stars in the solar vicinity. 
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It is known that intermediate mass stars (m < 8 M ) produce carbon and not 
oxygen. Carigi (1994) found that yields by Renzini & Voli (1981) can not explain the 
C/O increase with the metallicity in the solar neighborhood. Prantzos et al. (1994) 
suggested that C- yields higher that those predicted by Renzini & Voli could 
reproduce the observed rise in C/O. The other two sets of yields for LIMS dependent 
on Z that exist in the literature (van den Hoek & Groenewegen 1997, and Marigo et 
al. 1996, 1998) can not explain the recent rise of C/O in the solar vicinity. However, 
LIMS have an important role in the early evolution of C/O. 

6. CONCLUSIONS 

Based on basic chemical evolutionary models for the Galactic disk, assuming different 
sets of stellar yields dependent on metallicity, the following conclusions are reached: 

a) Different sets of stellar yields predict different C/O values. The recent rise of C/O 
with O/n is due mainly to massive stars. The early rise of C/O with O/H is due to 
both massive stars and LIMS. 

b) Models using stellar yields with stellar winds dependent on Z (Maeder 1992, 
Portinari et al. 1998) reproduce the rise of C/O with time shown by Hn regions, the 
Sun, and dwarf stars within the solar vicinity. 

c) Models considering Maeder 's yields are also successful in reproducing the C/O 
Galactic abundances and gradient determined from H II regions. 

d) The model assuming Padova yields reproduces C/O abundances but not C/O 
Galactic gradients from H II regions, because stellar winds oc Z - 5 become too strong 
for supermetallic stars. 

e) Models assuming yields from WW&WLW reproduce only the C/O abundances 
from B-stars, but fail to reproduce the other observational constraints. 

f) Since the C/H abundances in B-stars are lower than in disk dwarf stars of 
different ages in the solar neighborhood, while the O/H value is very similar, the 
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C/O value from these B-stars may not be a good observational constraint. 

g) The main difference between the sets of yields of massive stars arises because 
Geneva and Padova groups consider stellar winds dependent on metallicity while 
WW do not, and WLW only partially. One can then conclude that Z-dependent 
stellar winds must play an important role in the chemical enrichment history of the 
Galaxy. 

h) Observations within the solar neighborhood, the Galactic disk, as well as in spiral 
and irregular galaxies imply that C/O must increase with metallicity. Modeling 
including not only winds, but actually metal-dependent winds, is necessary to 
properly follow the chemical evolution of galaxies. 

i) To improve the agreement with the C/O Galactic abundances and the C/O 
evolution with metallicity, the present models call for a more complicated mass-loss 
rate law than Z - 5 , assumed by Maeder (1992) and by Portinari et al. (1998): 
M wind oc Z n , such that n > 0.5 if Z < Z Q , n ~ 0.5 when Z < Z & and 0.5 < n < 0.7 
for Z > Z Q . 

j) Models based on yields by Renzini & Voli (1981) predict less C and He than those 
by Marigo et al. (1996, 1998) and more than those by van den Hoek & Groenewegen 
(1997) 

k) The AY/AZ value decreases along the sequence of models MBCP - RV - HG. At 
the same time, it increases along the PCB - M92 - WW&WLW model sequence. 
1) The number of Hll regions with known C abundance is quite small, and the C/O 
gradient from H II regions might change with future C determinations in more H II 
regions. It is important to obtain C abundances from H II regions located at other 
galactocentric distances to determine the Galactic C/O gradient with higher 
accuracy. 

I would like to thank Manuel Peimbert for valuable comments and useful 
suggestions. I also wish to thank Laura Portinari for providing me the massive stars 
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yields of the Padova group. I am grateful to Don Garnett for supplying me the data 
for irregular and spiral galaxies. I also acknowledge the several excellent suggestions 
by the referee which have improved this final version. I thank Jesus Gonzalez for a 
thorough reading of the manuscript. 
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